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governed by intrinsic materials properties, and hence may be able to be reproducibly made. 
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I. INTRODUCTION 



Understanding the interface properties of strongly cor- 
related electron systems placed into inhomogeneous envi- 
ronments on the nanoscale combines the fields of strongly 
correlated electron systems and nanotechnology. The 
interface properties of strongly correlated systems will 
play an important role in determining the properties of 
devices made from these materials. One of the impor- 
tant interface properties is electronic charge reconstruc- 
tioni. At nearly all types of metal-semiconductor inter- 
faces, a so called Schottky barrier exhibits charge de- 
pletion in the doped semiconductor region close to the 
interface. Because the Fermi energy in the metal dif- 
fers from that in the semiconductor, mobile carriers in 
the semiconductor side of the barrier diffuse into the 
metal side until a static equilibrium is reached. Re- 
cent experiments conducted on strongly correlated ma- 
terials have shown a similar interface-induced charge re- 
construction in inhomogeneous nanostructures. Ohtomo 
and co-workers^ have fabricated atomically precise digi- 
tal heterostructures consisting of a controllable number of 
planes of LaTi03 (a correlated electron Mott-insulating 
material) separated by a controllable number of planes 
of SrTiOa (a more conventional band-insulating mate- 
rial). The experiment demonstrated an insulator- metal 
crossover near the interface due to electronic charge re- 
construction. The insulating heterostructure developed 
conducting channels near the interfaces for current par- 
allel to the planes. Okamoto and Millie 1 found, through 
a Hartree-Fock calculation, that the electronic charge re- 
construction leads to metallic behavior at the interface 
between the two insulators. This is because the mis- 
match of chemical potentials creates screened dipole in- 
terface charge reconstructions that are conducting due 
to the excess or deficit of charge. Varela and collabo- 
rators' 5 have provided evidence for extensive charge re- 
construction from manganite-cuprates heterostructures 
at a YBa2Cu307_a;/Lao.67Cao.33Mn03 interface, which 
also exhibits a metal-insulator crossover near the inter- 
face. From a theoretical standpoint, Nikolic, Freericks 
and Miller— developed a semiclassical approach to the 



Potthoff-Nolting 5 algorithm to create a dynamical mean- 
field theory description of electronic charge reconstruc- 
tion and applied it to superconductor-insulator-normal 
metal-insulator-superconductor Josephson junctions. 

In this work, we investigate the electronic charge re- 
construction of a doped Mott-insulator, where the recon- 
struction can dope parts of the system close to the insu- 
lating phase. Such a system may be realized in high tem- 
perature superconducting grain boundaries, where the 
electrically active grain boundary can lead to an elec- 
tronic charge reconstruction 6 . We take a semi-infinite 
ballistic-metal lead and couple it to another semi-infinite 
ballistic-metal lead through a strongly correlated barrier 
material of varying thickness. By adjusting the inter- 
action and the filling of the barrier material, we study 
interface properties for different types of materials. Here 
we emphasize the physics of the doped Mott-insulator. 

The organization of this paper is as follows: in Sec. 
II, we present a detailed derivation of the formalism and 
the numerical algorithms used to calculate the charge 
reconstruction of nanostructures. In Sec. Ill, we describe 
our numerical results. We end with our conclusions in 
Sec. IV. 

II. FORMALISM 

We apply the Potthoff-Nolting 5 approach to multilay- 
ered nanostructures, which involve translationally invari- 
ant x — y planes stacked in the longitudinal z-direction. 
We choose square lattice planes stacked along z direc- 
tion, so the lattice sites are identical to those of a simple 
cubic lattice. We have periodicity in the x and y direc- 
tions, but we allow inhomogeneity in the z-direction. All 
interactions are also translationally invariant within each 
plane, but can change from one plane to the next. The 
system is described by a spatially inhomogeneous many- 
body problem. Potthoff and Nolting described the idea 
of using a mixed basis for inhomogeneous DMFT: First, 
Fourier transform the x and y coordinates to wavevectors 
k x and k y but keep the z-component in real space. Then 
for each two-dimensional band energy, we have a quasi- 
one-dimensional problem to solve, which has a tridiago- 
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nal representation in real space, and can be solved with 
the so called quantum zipper algorithm 12 . 

Because of the translational invariance in each two- 
dimensional plane, we can describe the intraplane hop- 
ping via a two-dimensional bandstructure, which be- 
comes 



el (k x ,kj,) = -24[cos(k x ) + cosfkj,)], 



(1) 



for a square lattice plane, where t\ is the nearest- 
neighbor intraplane hopping on the erth plane and 
(k x , k y , 0) is the two-dimensional wave vector. 

For the interaction, we employ the Falicov-Kimball 
model£ which involves an interaction between conduction 
electrons with localized particles (/-electrons or charged 
ions) when the conduction electron hops onto a site oc- 
cupied by the localized particle. The Falicov-Kimball 
model has a non- Fermi liquid ground state in the metallic 
regime, because the electrons see static charge scatterers, 
which always produce a finite scattering lifetime, so there 
is no quasiparticle resonance, as seen in other strongly 
correlated models like the Hubbard ? - and periodic An- 
derson- model. It also has a Mott-type metal-insulator 
transition, that sets in when the correlation strength is 
large enough and the total number of the particles (local- 
ized plus delocalized) equals the number of lattice sites. 
This is because the energy cost for double occupation 
becomes too high if U is large enough, and the system 
becomes an insulator; the metal insulator transition can 
occur for systems without particle-hole symmetry, un- 
like the single-band Hubbard model which always has 
the Mott transition precisely at half filling-^. Our choice 
of using the Falicov-Kimball model is pragmatic, since 
the DMFT can be easily solved for this system and it de- 
scribes interesting metal-insulator transitions of strongly 
correlated materials. We expect the results in the insu- 
lating phase to resemble other correlated insulators, since 
the most important property of a correlated insulator is 
the size of its gap. On the metallic side, the Falicov- 
Kimball model is good for describing the crossover from 
ballistic to diffusive transport in dirty metals, but it is 
unable to describe the coherent quasiparticle formation, 
with a renormalized Fermi energy, seen in pure models 
like the Hubbard model. We feel it is nevertheless an 
interesting model to consider for examining systems near 
a Mott transition (especially since any experimental sys- 
tem will always have disorder, so the renormalized Fermi 
liquid will also disappear close to the Mott transition due 
to this disorder, and the actual metal-insulator transi- 
tion may be closer to the scenario of the Falicov-Kimball 
model). We consider spinless electrons here, but spin can 
be included by introducing a factor of 2 into some of the 
results; but note that it will modify the self-consistency 
relation for the Coulomb potential energy, which will not 
involve just a trivial change of the final converged results 
(the Coulomb potential energy will be doubled in magni- 
tude, but that doubled value goes into the local chemical 
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FIG. 1: Geometry taken for the classical electrostatics prob- 
lem. We show a blow up of two planes, a and a + 1. As- 
suming that the excess surface charge density on plane a is 
(p a — p^ lk )a = da and that the permittivity is e a (and sim- 
ilarly for the a + 1 plane), then the change in polarization 
at the interface between the two dielectric planes induces a 
polarization charge on the interface (denoted a po i) that leads 
to a discontinuous jump in the electric field halfway between 
the two lattice planes. Once all the fields are known, we in- 
tegrate the total field to get the electric potential. Note that 
the discontinuity in the electric field occurs at the midpoint 
between the two lattice planes. 



potential for the next iteration of the algorithm). In- 
troducing spin also changes the filling condition for the 
Mott-transition. 

The charge reconstruction at the interface leaves net 
charge on each plane. So the Hamiltonian needs to be al- 
tered to include this effect. We analyze this problem us- 
ing a semiclassical approach — the charge rearrangements 
are employed to determine classical electric fields, electric 
potentials and potential energies, which arc then input 
into the Hamiltonian, that is subsequently solved using 
quantum mechanics. We assume that the electric charge 
is uniformly spaced over the plane. Then the electric 
field is a constant, perpendicular to the plane. If the net 
charge density on plane a is p a — p h ^ lk and the permittiv- 
ity is e a , then the change in polarization at the interface 
between the two dielectric planes induces an additional 
polarization charge on the interface that leads to a dis- 
continuous jump in the electric field halfway between the 
two lattice planes. The local field created by the crth 
plane has magnitude 



\E\ = 



ea(p 



n bulk\ 



2e a 



(2) 



in which a is the lattice spacing along z direction. Once 
the total fields at each plane are known, we integrate 
them to obtain the electric potentials. The discontinuity 
in the electric field occurs at the midpoint between the 
two lattice planes when there is a change in permittivity. 

Performing the integration yields the potential energy 
due to the Coulomb interaction of the electronic charge 
reconstruction as 
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, J2y= a +Ai e Schoth) + eschoth ~ 1)], P>a 
V = -J2(p a -Pa lk )^<{o, a & = (3) 

E°=(j+i \[eSchot{l) + eschoth - 1)], P < a 



where we define the symbol es c hot{ct) = e 2 a/2e a , which 
is related to the screening length in a particular medium. 
This parameter has the units of an energy multiplied by 
an area; the product of eschot with the local DOS has 
the units of the inverse of a length, and this is what 
determines the decay length of the charge profile. In our 
calculations we will set eschot to a particular value to fix 
the screening length. 

The Coulomb potential energies are input into the 
Hamiltonian as follows: they are treated by shifting the 
chemical potential /i — ► ji— V a on each plane, depending 
on what the Coulomb potential energy is. The Falicov- 
Kimball Hamiltonian with electronic charge reconstruc- 
tion is then 

aiflj oti 

- ^(^ + AS Fct -V4)cLc m , (4) 

oti 

in which a and (3 represent different planes and i and 
j represent lattice sites on the respective planes. The 
chemical potential fi is fixed by the bulk value of the 
leads. AEp a is the mismatch of the center of the bands 
between the leads and the barrier (AEp a = in the 
leads). By adjusting the value of AEp a , we control the 
relative shift of the bands inside the multilayered nanos- 
tructure; indeed, the electronic charge reconstruction oc- 
curs because AEf is adjusted to describe the chemical 
potential mismatch between the two materials at some 
given temperature. Note that AEf is a fixed constant 
that does not change with the temperature. 

To analyze the many-body problem, we use a Green's 
function technique. The imaginary-time Green's func- 
tion, in real space, is defined by 

G m p 3 {r) = -(T T c m (r)4(0)). (5) 

for imaginary time r. The notation (O) denotes 



TrT T exp(-/3[H - [iN])0/Z, where Z is the partition 
function Z = Trcxp(— f3[H — ^N]). The operators 
are expressed in the Heisenberg representation where 
0(t) = exp(r[H ~ nN})0 exp(-r[H - nN}). The symbol 
T T denotes time ordering of the operators, with earlier r 
values appearing to the right and (3 is the inverse tem- 
perature, /3 = l/T. We will work with the Matsubara 
frequency Green's functions, defined for imaginary fre- 
quencies iu> n = iiiT(2n +1). These are determined by a 
Fourier transformation 

G a i0j{iu) n ) = / dTe lUnT G aiP] {T). (6) 
Jo 

We can write the equation of motion for the Green's func- 
tion in real space, which satisfies 

G aif3j( i(J n) = {iw n + fl + AE Fa - V a )8 m pj 

— T, a (i(jJ n )5 a ij3j + t a ipj. (7) 

Now we go to a mixed basis, by Fourier transforming in 
the x and y directions, to find 

G~p(k, iuj n ) = [iu> n + (m + AE Fa - V a ) - T, a (icj n ) (8) 

with Y, a {iuj n ) the local self-energy for plane a and k the 
two-dimensional planar momentum. Finally, we use the 
identity ^ 7 G Q7 (k)G~^(k) = S a p to arrive at the start- 
ing point for the recursive solution to the problem, 

S a p = G a p(k,iu; n ) 

x [iu n + + AE F p - Vp) - T,p(iu n ) - c J(k)] 
+ G a p-i(k,iuj n )tp-ip + G a p +1 (k,iuj n )tp +1 p. (9) 

It turns out that there is a straightforward procedure to 
determine G a p from this equation of motion. It is called 
the quantum zipper algorithm. We start with (3 = a, 
which can be used to find the Green's function via 



G aa (k, iu n ) = fi r< 7T~ n F TT - n ■ (10) 

iuj n + (p + AEp a - V a ) - X a (iw n ) - ei(k) + %^ig^l ta _ lQ + ^gffi k+la 

Next, we consider the equations with which can be put into the form 

Gota— m+1 (k, 2(J n )£ a _ m+la— m ■ ./.Am \ v / * \ II /i \ , ^ocol— rn — 1 ( ^ , i U) n )t a — m _ i a _ rn 

^ 7j~ . , = iu„Hv+ AE F*- m -V a - m )-Y, a _ m (iu n )-e a _ m {lt)+ — , 

(11) 
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for m > 0, with a similar result for the recurrence with 
to < 0. In these equations, we have used E Q to denote 
the local self-energy E QQ on plane a. We define the left 
function 



Gcta — m+l (k, iw n )f Q - m ^i Q -f 



(12) 



and then determine the recurrence relation from Eq. (|llp 



L a - m (k, iuj n ) = iu) n + ((*+ AE Fa _ m - V a - m ) ~ 

t a— mot— m— 1 « a— m— Ice— m 



^ot—m 



L (k, iw n ) 



•(13) 



We solve the recurrence relation by starting with the re- 
sult for L_ oc , and then iterating Eq. [)13p up to to = 0. 
Since we must have a finite set of equations for an ac- 
tual calculation, we assume we have semi-infinite metallic 
leads, hence we can determine L_oo by substituting L_oo 
into both the left and right hand sides of Eq. (fT"3|) with 
AEpa + V a =0, which produces a quadratic equation 
for L_oo that is solved by 



L_ oc (k,iw„) 



iu n + (i - S_oo(iw„) - eloJk) 



'[iw„ + M - S_co(iw„) - eloJk)] 2 - 44=00 . ^ 

± o -( 14 ) 



This determines the left functions far from the interface. 
The sign in Eq. (|14p is chosen to yield an imaginary part 
less than zero for iio n lying in the upper half plane, and 
vice versa for iw n lying in the lower half plane. When we 
are sufficiently far from the interface, the Green's func- 
tions will be essentially the same as the bulk and hence 
the L a functions will equal £-oo- In our calculations we 
allow L a to differ from £_oo only for the thirty planes 
closest to the interface on either side of the barrier. This 
means we start the recurrence relation in Eq. (| 1 3[) with 
ol — in — 1 being the thirty first plane to the left. Then 
all subsequent L^s are allowed to vary until we reach 31 
planes to the right of the barrier, where we assume L a 



J 



becomes a constant again. This approach is accurate, 
when the system heals to its bulk values within those 
thirty planes on either side of the interface. If this heal- 
ing has not occurred, then one needs to include more 
planes before one terminates the problem with the semi- 
infinite bulk solution. 

In a similar fashion, we define a right function and a 
recurrence relation to the right, with the right function 
being 



-Ra+m(k, lUJ n 



Gaa+m— 1 (k, ^^n)^a+m— la+r. 



(k, iu n ) 

and the recurrence relation satisfying 
# a+m (k, iuJ n ) = iu n + (^ + AE Fa+m - V a ^ 



(15) 



E Q + m (iw n ) 



e 11 

ct-f-m 



(k)- 



a+rn+l^a+m+la+m 



(16) 



-Ra+m+l(k, iljJ n ) 

We solve the right recurrence relation by starting with 
the result for Roo , and then iterating Eq. (|16p up to to = 
0. As before, we determine Roc by substituting Roo into 
both the left and right hand sides of Eq. (|16p . where 
AE Fa -V a = 0, 



iuj n + £oo(icj„) - eio(k) 



± 



\j [iw n + - £oo(«^«) - elo(k)] 2 - At 2 ^ 



(17) 



The sign in Eq. (fT7|) is chosen the same way as in Eq. (TTj 
In our calculations, we also assume that the right func- 
tion is equal to the value Roo found in the bulk, until 
we are within thirty planes of the first interface. Then 
we allow those thirty planes to be self-consistently de- 
termined with R a possibly changing, and we include a 
similar thirty planes on the left hand side of the last in- 
terface, terminating with the bulk result to the left as 
well. Using the left and right functions, we finally obtain 
the Green's function from 



G aa (k,iLu n ) = 



1 



L a (k, iu> n ) + R a (k,iw n ) - [iu„ + + AE Fa - V a ) - S Q (iw n ) - e«(k)] 

I 



(18) 



where we used Eq. ([T3J) and Eq. (jl6|) in Eq. JTO]). This 
technique for determining the mixed-basis Green's func- 
tion is called the quantum zipper algorithm^; it has been 
modified here to treat the electronic charge reconstruc- 
tion. 

The local Green's function on each plane is then found 
by summing over the two-dimensional momenta, which 
can be replaced by an integral over the two-dimensional 



density of states (DOS) since all the momentum depen- 
dence in the algorithm is in terms of ei: 



G aa (iw n ) = / delp 2d {el)G aa (el7Lu n ), (19) 
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with 



2d( II \ 



\—K I 1 - . 
2-kHW \ \ 



1 



(el? 
(4ti) 



(20) 



and K(x) the complete elliptic integral of the first kind. 

The DMFT algorithm starts with a self-energy on each 
plane, which is usually chosen to be zero. Next, we use 
the quantum zipper algorithm to find the local Green's 
function on each plane. This step is the inhomogeneous 
nanostructure equivalent to the Hilbert transform, which 
is used in bulk DMFT. Once the local Green's function 
is known on each plane, we extract the local effective 
medium via 



ensure stability of the solutions. If the potentials are up- 
dated too quickly, then the system of equations migrates 
to an unphysical fixed point, or limit cycle. Typically, 
we use 0.99 of the old potential (or more) and 0.01 (or 
less) of the new potential in each averaging step. It usu- 
ally takes more than a thousand iterations to reach the 
desired convergence, but because our impurity solver is 
for the Falicov-Kimball model, the calculations can still 
be completed rapidly. 



III. NUMERICAL RESULTS 

Previous work on the multilayered nanostructure at 
half filling has shown interesting properties^. For exam- 

bulk\ 



GodiuJny 1 = G^iLUn)- 1 + T, a {iuj n ), (21) pie, if the barrier is metallic, the net charges {p a - p™ (fc ) 



for each plane. Next, we need to solve the local impurity 
problem for the given Hamiltonian on the ath plane with 
the given effective medium: 



G a {iu n ) = (1 - wi)G 0a (iu n ) + 



Wl 



G a{iUn) 1 - U a 



> (22) 



with wi the /-electron density. This will produce a new 
local Green's function for each plane, and then a new 
self-energy for each plane via Eq. pTjl . 

The electronic charge on each plane is calculated by 
summing the Green's functions over all Fermionic Mat- 
subara frequencies on the imaginary axis, 



3 Q = - + T^G a (iuj n ). 



(23) 



Since G a (iuj n ) behaves like l/iu) n at large n, we regu- 
larize the summation by subtracting T^2 n l/[iu>„ — p + 
ReS Q (iti; nmax )] from G a (iui n ) to speed up the conver- 
gence. Re£ Q (iw„ max ) is the real part of the self energy 
at the highest Matsubara frequency (n max ) we use in 
our calculation. Typically n max satisfies 2-7rTn max = 
30. Since T^-^b^ = ^nh{/% 



iuj n — z^+RcSq, (it^Ti max ) 

ReE a (iu; nniax )]/2}/2, Eq. |g3J) becomes 



pa 



G a (iLu n ) 



1 



+ — tanh . 



iuj n — p + RbEq, {w nms 
I (3[p — ReS ct (iti; rimax 



(24) 



Note that our regularization scheme requires that 
TJ2 n G a (iuj n ) = (4 (O)ca(O) - c Q (0)4(0)>/2. Once the 
change in the charge density is known, we can calcu- 
late the electrical potential V a through Eq. (J3J). This is 
then added to the chemical potential to determine the 
electrochemical potential at each plane. Then we iter- 
ate the DMFT algorithm with the new electrochemical 
potential until the self-energy at every plane converges 
and the potentials no longer change. The iterative equa- 
tions need to be solved with an averaging procedure to 



on every plane divided by the amount of the band shift 
(AEf) exhibit scaling: the curves of the charge recon- 
struction with different band shifts essentially lay on top 
of each other with deviations only occurring very close to 
the interface. While this effect is most likely due to the 
flatness of the simple cubic lattice density of states near 
the band center, the result was still rather striking in the 
data. 

In this contribution, our focus is on the doped Mott- 
insulating phase. The doped Mott-insulator-metal inter- 
face is in many ways similar to the Schottky barrier. The 
barriers are both insulating; the doped carriers move to 
the interfaces as the result of the band mismatch, causing 
depletion of carriers close to the interface inside the bar- 
riers. However, the doped Mott-insulator is of particular 
interest because the band gap and the small amount of 
doped charges generate interesting physics that may be 
useful for designing future devices^. 

We reduce the number of parameters in our calcula- 
tions by assuming all of the hopping matrix elements arc 
equal to t for nearest neighbors. While not required, it al- 
lows us to reduce the number of parameters that we vary 
in our calculations so that we can focus on the physical 
properties with fewer calculations. The hopping scale t is 
used as our energy scale. As described above, we include 
30 self consistent planes in the metallic leads to the left 
and to the right of our barrier, which is varied between 15 
and 30 planes in our calculations. The screening length, 
as determined by the parameter eschot( a )i is about 2.2 
lattice spacings. We choose eschot(a) — 0.4 through- 
out the device (in both the leads and the barrier). The 
band shift AEf can be divided into two parts, Ef(T) 
and Pbarrier(T)- I n the following analysis, we plot our 
result with respect to E f {T), which is AE f -p} b ^ k rier (T). 
The reason for doing this is that if we shift the band of 
the barrier to the amount of the bulk chemical potential 
of the barrier, the chemical potentials of the lead and 
the barrier should be automatically aligned with each 
other. Shifting the band less than that (Ef < 0), nega- 
tive charges should be attracted to the barrier in order 
to reach equilibrium; more than that {Ef > 0), positive 
charges should be attracted to the barrier. (Note that 
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FIG. 2: Density of states of the central plane in the bar- 
rier versus frequency for U = 12, wi — 0.75, Ef = and 
p e = 0.251. The temperature is T = 0.25. With strong in- 
teraction, the density of states has a large gap between the 
bands. Note the size of the band gap is roughly 5.0 and if 
the carriers are drawn from the band across the gap, \Ef\ 
has to be larger than the size of the band gap. This DOS 
is indistinguishable from the bulk DOS on this linear-scale 
plot. The density of states does not change with T, except 
for the location of the zero on the frequency axis, due to the 
temperature dependence of the chemical potential (as T — » 0, 
the origin lies just slightly above the lower band edge of the 
upper Hubbard band). 
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because the local Green's functions change near the in- 
terface, even if Ef(T) = 0, there will always be a small 
charge reconstruction in the general case.) At half filling, 
the curves are particle-hole symmetric between positive 
and negative Ef. At half filling, the bulk chemical po- 
tential is independent of the temperature, so changing 
the temperature will not have effects on the mismatch of 
the bands. At fillings other than half filling, the chemi- 
cal potential is dependent on temperature. So changing 
temperature can have a similar effect as changing Ef. 
When the temperature is low, the effect of different tem- 
peratures on the charge reconstruction is small because 
of the large band gap and also because the change in 
the bulk chemical potential is small at low temperatures. 
In our calculation, we set T = 0.25. If we take a reason- 
able energy scale for our system, such as a noninteracting 
bandwidth of 3 eV, then t = 0.25 eV, and the tempera- 
ture corresponds to 750 K. 

We tune the Falicov-Kimball interaction in the barrier 
planes to U — 12, which lies well into the Mott-insulating 
regime. We set the conduction electron filling at 0.251 
and the density of the /-electron w\ at 0.75. Hence, we 
are considering a slightly doped particle-hole asymmet- 
ric Mott insulator. The bulk DOS of the barrier planes 
is shown in Fig. [2] The integration of the DOS over 
frequency in the lower band is 0.25 and the integration 



FIG. 3: Charge deviation in a multilayered nanostructure 
with a barrier thickness of 20. The barrier is slightly doped 
with electrons so that the conduction electron filling is 0.251 
[/Cner(r = 0.25) = 7.046]. The different lines denote dif- 
ferent Ef values, as detailed in the legend. Note the charge 
deviation pattern with (a) negative Ef is different from the 
pattern with (b) positive Ef. With positive Ef, the magni- 
tude of the charge deviation on the interface increases with 
the magnitude of the shift. With negative Ef, if the mag- 
nitude of the shifts are small, the magnitudes of the charge 
deviation do not change too much due to the presence of the 
Mott gap. (Color on-line.) 

of the DOS in the upper band is 0.75. Because of the 
strong interaction, there is a gap with a size of roughly 5 
between the lower band and the upper band. 

The existence of a large band gap in the strongly cor- 
related doped insulator dramatically changes the charge 
reconstruction from the metallic case. The device ex- 
hibits asymmetry with regard to positive and negative 
band shifts as shown in Fig. [31 With positive band shifts, 
the carriers can easily move from the leads to the barrier 
until the electric field generated by electronic charge re- 
construction strikes a balance with the mismatch of the 
bands. Larger band shifts cause larger electronic charge 
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FIG. 4: Electric potential (V a ) as a function of Ef. The 
conduction electron filling is 0.251 and the interaction is U = 
12. Note the drastic increase of the potential at the center of 
the barrier from E F (T = 0.25) = -4 to E F (T = 0.25) = -6 
(recall the bandgap is equal to 5). (Color on-line.) 



reconstruction at the interface. The situation is more 
complicated for negative band shifts. 

If the shift is small, the doped carriers are completely 
drained from the barrier (p e — > 0.25) but the band gap in 
the Mott insulator stops further diffusion of the carriers 
into the leads. The charge profile in this case is similar to 
that of a Schottky barrier. In the n-type Schottky bar- 
rier, electrons in an n-type semiconductor can lower their 
energy by traversing the junction. As the electrons leave 
the semiconductor, a positive charge, due to the ionized 
donor atoms, stays behind. This charge creates a nega- 
tive field and lowers the band edges of the semiconductor. 
Electrons flow into the metal until equilibrium is reached 
between the diffusion of electrons from the semiconduc- 
tor into the metal and the drift of electrons caused by the 
field created by the ionized impurity atoms. The differ- 
ence is that in a Schottky barrier, a region in the semi- 
conductor close to the junction is depleted of mobile car- 
riers while in our multilayer nanostructure, the carriers 
are completely drained from the doped barrier because 
of the small barrier thickness. (Note in Fig. [3ja), with 
small shifts, the charge deviations in the barrier planes 
uniformly have the value of —0.01 except for the planes 
close to the interfaces, which means all the doped carriers 
move to the interfaces.) When the shifts get larger, the 
carriers close to the interface begin to move to the leads. 
Thus holes are created near the interface in the barrier 
side. When the band shift gets large enough, the band 
gap can no longer stop the diffusion of the carriers from 
the barrier to the leads. The charge deviation curve at 
large negative band shifts bears a similar shape to the 
curve at positive band shifts, although the maximum net 
charges on the interface are significantly less due to the 
effect of the Mott-insulating band gap. 

To better understand the effect of the band gap, Fig.0] 
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FIG. 5: Charge deviation in a multilayered nanostructure 
with a barrier thickness of 20. The barrier is slightly doped 
with holes so that the conduction electron filling is 0.249 
[tJ-balrieriT = 0.25) = 2.488]. The different lines denote differ- 
ent Ef values, as detailed in the legend. The charge deviation 
pattern in this case is the opposite of that with the 0.251 fill- 
ing. With negative Ef [panel(b)], the magnitude of the charge 
deviation on the interface changes with the magnitude of the 
shift. With positive Ef [panel (a)], if the magnitude of the 
shift is small, the magnitude of the charge deviation does not 
change much. (Color on-line.) 



shows the electric potentials generated by the charge re- 
construction with Ep = —4.0 and Ep — —6.0 respec- 
tively. The nearly five fold increase of the potential at 
the center of the barrier is because with Ep = —6.0, the 
band gap can no longer hold the carriers in the lower band 
from moving to the leads, resulting in more net charges 
along the interface. This is quite different from the case 
where the barrier is metallic, where an increase of Ef 
generically results in a smooth increase of the potential. 

We find opposite results if the carriers are holes in- 
stead of electrons. Fig. [5] shows the electronic charge 
reconstruction for the doped Mott insulator with a con- 
duction electron charge density p e — 0.249. In this case, 
negative band shifts attract holes from leads to the bar- 
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FIG. 6: Charge deviation as a function of interaction U. The 
conduction electron filling is 0.251 and the shift Ef is fixed 
at -2.0 [Mbarrier(r = 0.25) = 1.792 for [7 = 6, ^ h b l l r k rier {T = 
0.25) = 3.250 for U = 8 and y%£ T k rier {T = 0.25) = 5.121 for 
U = 10]. Note the magnitude of the charge deviation on the 
interface drops as U increases. However if U is large enough, 
further increase of U will not result in a further decrease of 
the magnitude of the charge deviation. (Color on-line.) 



rier. Larger shifts create more net charges along the in- 
terface. Small positive band shifts can only drain holes 
from the lower band of the barrier due to the band gap. 
For large positive band shifts, the charge deviation curve 
returns to the shape of those with negative band shifts. 

Next we examine the electronic charge reconstruction 
for a given amount of band shift with different strengths 
of the interaction. We choose electrons as the carriers 
Oe = 0.251). The band shift E f is fixed at -2.0. Fig. [6] 
shows the different charge profiles as U changes from 6 to 
10. For small U, the gap of the Mott insulator is small. 
The electrons inside the barrier can easily diffuse into the 
leads. We see a large electronic charge reconstruction at 
the interface. As U gets larger, it becomes harder for 
the electrons to move into the leads due to the larger 
band gap. The net amount of the charge deviation at the 
interface becomes less. With strong enough interaction, 
only the doped carriers can be drawn from the barrier. 
The curve then resembles the one with [7 = 12. 

In our model, the screening length eschot can be 
changed in the leads as well as in the barrier. We tried 
several different values of eschot ranging from 0.1 to 2.0. 
With a larger screening length, the magnitude of the 
charge reconstruction increases, essentially because the 
electric field generated by the charge reconstruction is 
stronger. There are more net charges on the planes close 
to the interface. Planes deeper inside the leads are less 
affected by the change in the screening length. The ef- 
fect of eschot is less obvious in the doped Mott insulator 
case with large interaction. The total number of carriers 
is limited in this case. A longer screening length can- 
not move more carriers into the leads, because they are 



depleted near the interface and the field is not strong 
enough to start moving charge out of the lower Hubbard 
band. 

The thickness of the barrier is not playing an impor- 
tant role in these multi-layered structures (for moderately 
thick barriers). We adjusted the thickness of the barrier 
from 10 to 20 planes. More planes in the barrier obvi- 
ously provides more carriers, generating more net charges 
on the interface. But the properties of the central part 
of the barrier remain similar. 



IV. CONCLUSION 

In this paper, we applied a generalized DMFT to inho- 
mogeneous systems to calculate the self-consistent many- 
body solutions for multilayered nanostructures with bar- 
riers that can be adjusted to go through the Mott transi- 
tion. We developed the computational formalism based 
on the algorithm of Potthoff and Nolting and the Falicov- 
Kimball model that can calculate the charge reconstruc- 
tion of the multilayered nanostructures with variable bar- 
rier fillings. We focused our study on the doped Mott 
insulator as the barrier material. We found interesting 
results that came out of this analysis. 

First of all, the scaling effect (the charge reconstruc- 
tion divided by the value of Ef remain roughly the same 
across the system) at half filling in the metal phase is 
no longer valid at fillings other than 0.5. Second, the 
symmetry of the charge reconstruction with positive and 
negative Ef hi the metal phase is broken in the Mott in- 
sulating phase. With electrons as carriers in the doped 
barrier, small negative Ef can only draw the doped car- 
rier to the leads. If the value of \Ef\ is smaller than the 
size of the band gap of the Mott insulator, the net charges 
on the interface do not increase much with the increase of 
\Ef\ and the central part of the barrier remains drained 
of carriers. If the the value of \Ef\ rises above the size 
of the band gap of the Mott insulator, electrons from the 
lower band start to move to the leads. The charge profile 
changes drastically. With holes as carriers in the doped 
barrier, the positive region of Ef is affected by the band 
gap. So the charge reconstruction is opposite to the elec- 
tron doped barrier with respect to Ef. Third, the charge 
reconstruction of the doped Mott insulator case does not 
change much with changing parameters like eschot or the 
thickness of the barrier. The central planes of the barri- 
ers are always drained of carriers. And the amount of the 
net charges on the interface is largely determined by the 
number of carriers inside the barrier with the strength of 
the interaction smaller than the band gap. That means 
physical properties of devices made in this structure may 
be stable even when the parameters of the barrier mate- 
rial mentioned above have some deviation. So the next 
step is to calculate the transport properties or the resis- 
tance of this multilayered nanostructure with the doped 
Mott insulating barriers and to compare the results from 
different parameter sets. This requires generalizing the 
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codes to the real axis, which is beyond the scope of this 
work. With small band shifts, we anticipate that the 
resistance should not change much because the charge 
reconstruction is dominated by the amount of carriers 
in the barrier. Because the carriers in the barrier are 



completely drained, the resistance should be large. With 
a large band shift, the resistance should be significantly 
lower since charges from the lower band now move to the 
lead, and the whole structure is more conductive as a 
result (even though the Coulomb potentials are larger). 
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